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Abstract 

In this Reply to the preceding Comment of Klimchitskaya and 
Mostepanenko (cf. quant-ph/0703214), we summarize and maintain 
our position that the Drude dispersion relation when inserted in the 
Lifshitz formula gives a thermodynamically satisfactory description 
of the Casimir force, also in the limiting case when the relaxation 
frequency goes to zero (perfect crystals). 
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The main concern of the authors of Ref. pQ is that in their opinion the 
Drude dielectric function is not valid if the metal is a perfect crystal, with 
no impurities at all. By contrast, if the metal is imperfect implying some 
relaxation, the authors seem to agree with our standpoint as expressed in 
Ref. pj: The Drude dispersion relation is adequate to describe the medium 
and there is no conflict with the Nernst theorem in thermodynamics. The 
Drude relation implies that there is a finite relaxation frequency, which we 
call v. We find it satisfactory that now there seems to be full agreement 
as far as imperfect crystals are concerned. The problem remains for perfect 
crystals only. 

We can agree with the claim of Ref. pQ that perfect crystals is an ideal- 
ization that forms much of the basis of solid state physics. From a physical 
viewpoint taking the limit u ^ is however in the Casimir context a delicate 
and not unique operation. There are several ways to argue that our results 
found via use of the Drude relation are valid in general. 

1 In practice, all solids have some relaxation. This follows also from the 
Kramers-Kronig relations, for any material having frequency disper- 
sion. Ideal crystals can be seen as the limit of no relaxation v — > 0. 
We have shown very accurately, both analytically and numerically [2] 
- cf. also Refs.[3l H] - that for v > the entropy S = at T = 0. From 
this we draw the conclusion that in the limit v — > the value S = is 
preserved. But the entropy decreases more and more rapidly when T 
increases from T = 0, and this decrease becomes a discontinuous jump 
in the limit v — > 0. The variation of the entropy as a function of T for 
two copper half-spaces is nicely illustrated in Fig. 3 in the QFEXT03 
Proceedings article of Sernelius and Bostrom [5], for a finite value of v. 

2 Further, merely taking the limit v —>■ is not fully satisfactory from 
a physical viewpoint. It should be noted that on the molecular level 
the permittivity becomes non-local, e = e(u,k), where k is the mag- 
nitude of the wave vector k. Then e becomes finite for finite k. Only 
the special case e(0, 0) is infinite, and this ("measure zero") can not 
be expected to give a finite contribution to the Casimir force. In inde- 
pendent studies Svetovoy and Esquivel [6] and Sernelius [7] find that 
the TE zero mode should not contribute when spatial dispersion is 
carefully taken into account. Because of nonlocal effects, the former 
of these state, "the question does [v\ go to zero or have some residual 
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value at T — > becomes unimportant" and they find that the TE zero 
frequency contribution must indeed be zero to satisfy Nernst's theorem. 

3 A TE zero frequency mode is not a solution of Maxwell's equations and 
should not occur for that reason. This point is discussed in more detail 
in Refs. [I] and [5]. (See especially Sec. Ill of Ref. [8].) 

4 Introducing a TE zero frequency mode for ideal crystals would imply 
that a medium with v = would behave differently from a real medium 
when taking the limit v — > 0. In the former case with the TE-mode 
included the entropy increases slowly from S — when T increases, 
while in the latter case it decreases rapidly to negative values for a 
real medium with v > 0. Thus in the limit v — > the result of the 
latter is different from the former. Due to this the use of the latter for 
v — > and the former for v = would create a discontinuity in physical 
properties at v = that we find unphysical. 

In connection with Eq. (3) in Ref. [1] it is argued that v(T) <C Cm CO — 
27rkmT/h (m = 1,2,3, •••) holds as temperature decreases. Surely, 
this is so for a perfect lattice in which relaxation follows the Bloch- 
Griineisen formula all the way to zero temperature. We consider real 
metals with impurities, however, in which v reaches some residual value 
at low temperatures. At temperatures so low that relaxation is dom- 
inated by impurities, their equation (3) no longer holds, but for suffi- 
ciently small T, Eq. (1) does. The authors of [T] explain in great detail 
why, for this reason, our calculated low temperature correction is nu- 
merically inaccurate. This is true, but immaterial, since our analysis is 
phenomenological and does not aspire to provide exact numerical esti- 
mates. The numerical calculation presented was motivated by a wish to 
confirm the correctness of our analytical result, which it convincingly 
does. 

The question becomes how to define the perfect lattice at zero temper- 
ature since the limits i/ — > and T — > do not commute. One may 
put v — right at the outset, or take the limit of vanishing relaxation 
frequency subsequent to calculating the free energy. For reasons dis- 
cussed above we suggest the latter procedure is the more appropriate. 
If equation (3) of pQ is to hold at all temperatures, the former proce- 
dure must be employed, which results in a different prediction of the 
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Casimir force and fulfillment of Nernst's theorem by inclusion of the 
TE zero frequency mode. 

Let us sum up: The basic assumption underlying our calculation in 
Ref. [2] was that the relaxation frequency v stays finite at any temper- 
ature including T = 0. At low temperatures we expect v to be smaller 
than at room temperature; this changes our results quantitatively but 
not qualitatively. On the basis of this assumption, we showed that 
there is no thermodynamic inconsistency associated with the Drude re- 
lation. Now, simply setting v = makes our formalism to some extent 
indeterminate. In particular, this indetermination persists if v = 
for T = 0. Physical arguments, like the requirement of continuity of 
behavior with respect to different values of input parameters, suggest 
strongly, however, that the results obtained from the Drude relation 
are quite general. 

References 

[1] G. L. Klimchitskaya and V. M. Mostepanenko, preceding Comment, 
Phys. Rev. E (EFK1007). 

[2] J. S. H0ye, I. Brevik, S. A. Ellingsen, and J. B. Aarseth, Phys. Rev. E 
75, 051127 (2007). 

[3] I. Brevik, S. A. Ellingsen, J. S. H0ye, and K. A. Milton, Preprint arXiv: 
0710.4882. 

[4] J. S. H0ye, I. Brevik, J. B. Aarseth, and K. A. Milton, Phys. Rev. E 67, 
056116 (2003). 

[5] Bo E. Sernelius and M. Bostrom, in Quantum Field Theory under the 
Influence of External Conditions", edited by K. A. Milton (Rinton Press, 
Princeton, NJ, 2004), p. 82. 

[6] V.B. Svetovoy and R. Esquivel, Phys. Rev. E 72, 036113 (2005) 

[7] Bo E. Sernelius, Phys. Rev. B 71, 235114 (2005). 

[8] J. S. H0ye, I. Brevik, and J. B. Aarseth, Phys. Rev. E 63, 051101 (2001). 



4 



